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Abstract 

Using the recent CLEO measurement of Br{b 57), we find that the branching 
ratio of b ^ sg cannot be larger than 10% in two Higgs doublet models. The 
small experimental value of Br{b evX) can no longer be explained by charged 
Higgs boson effects. 
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It is well known that the process 6 ^ 57 is extremely sensitive to new physics beyond the 
standard model, in particular, that containing a charged Higgs boson. In 1993, the CLEO 
collaboration placed an upper limit of Br{b —>■ 57) < 5.4 x 10^"^ on the inclusive branching 
ratio 1^ , which has inspired a large number of studies of this decay in various models for new 
physics [0 . Stringent constraints are obtained. Recently, CLEO has measured the inclusive 
branching ratio to be g] (2.32 ± 0.51 ± 0.29 ± 0.32) x 10-^ in short, (2.3 ± 0.7) x 10-^ 
corresponding to the 95% confidence level range of 

1 X 10"^ < Br{b ^ S7) < 4 X 10"^ (1) 

In the standard model, there is a surprisingly large QCD enhancement of 6 ^ S7 Q 
amplitude. This has stimulated intense efforts in calculating QCD corrections to leading 
order (LO), as well as partial calculations to next-to-leading order (NLO) 0,^. The CLEO 
result of eq. (1) is not far from the prediction of Br{b — > S7) in the standard model. This 
implies that not much room is left for new physics contributions to other b ^ s transitions 
such as 6 ^ sg, where the emitted gluon is "on-shell". 

It was shown in ref. |^ that for some choices of parameters in two Higgs doublet models, 
charged Higgs boson effects may enhance the decay branching ratio of b ^ sg beyond the 
10% level. Grz§dkowski and Hou |§] have pointed out that iib ^ sg rate is at the (10 — 20)% 
level, the discrepancy on Br{b —>■ evX) between experimental measurement (10.7 ± .5%) 
and theoretical expectations (> 12% in the standard model) |10| could be resolved. It is 
therefore of interest to check whether the possibility of Br{b sg) ~ 10% still holds once 
one includes the constraint imposed by eq. (1). 

In this report we focus on two Higgs doublet models (2HDM). These models are very 
simple extensions of the standard model, yet they exhibit some of the characteristics of 
a more complicated scalar structure typical of most theories beyond the standard model. 
We will consider the two distinct models (I and II) that naturally avoid tree-level FCNCs 
|n|. In Model I, one doublet (^2) couples to all fermions and the other (0i) decouples 



from the fermion sector. In Model II, 02 couples to up-type quarks while 0i couples to 



down-type quarks. This type of model occurs in minimal realization of supersymmetry, or 

in models with a Peccei-Quinn symmetry |T^. The major non-standard feature of these 

models is the appearance of extra physical scalar fields. We consider only the effect of 

charged Higgs bosons. Two parameters are sufficient to account for the additional effects. 

We take these to be rriH, the mass of the charged Higgs boson, and ^ = "^1/^2, the ratio of 

the vacuum expectation values of 0i and 02- Note that ^ = 1/ tan/?, as is commonly used 

in super symmetric models. In 2HDM, the charged Higgs couple to quarks with the same 

quark mixing matrix as the standard charged current. 

The standard model calculation for h — >■ S7, including up to date QCD corrections, can 

be found, for example, in Refs. |^,|^. The inclusive branching ratio is given by 

Y{h — * 57) 

Brih 57) = — — Brih —>■ cev) 

T[b ceu) 

_ \V:M^ 6a 1 ic^e//(^)|2 ^^(^ ^ (2) 



\Vch\^ 7rf{mc/mb) Q.{mt / Mw , fJ^) 
where the Wilson coefficient 



C!;"{fi) = T^^CriMw) + - (r/^ - r?i) Cs{Mw) + C2{Mw) <^^v'% (3) 



i=l 



includes short distance effects at Mw scale, while perturbative QCD effects are accumulated 
when running down to the physical scale /i, with 77 = as{Mw)/<ys{f^)- In eq. (0), the phase 
space factor f{z) is given by 

f{z) = l- + 8/ -z^ - 24z^ In z, (4) 

and the quantity Q{z) contains the 0{as) QCD corrections to the semileptonic decay rate 



T3| , p^ and is given by 

(5) 



The scheme-independent numbers and bi are given by [|15 



a, = (fill, -fff|, -f, -3^, -0.6494, -0.0380, -0.0186, -0.0057), 
^^=(i' i' i' -i' 0.4086, -0.4230, -0.8994, 0.1456), (6) 
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respectively. Defining x = mf/M^^, h = mf/Mfj, the coefficients Ci{Mw) are 



C^iMwY = 1, 



A{x)+e 



B{h) - -A{h) 
E{h) - ^Dih) 



for Model I, and 0,|T6| 



c2{Mwy' 



C-riMwY' = --Aix) 



-Dix 



-B{h)-\eA{h), 



-E{h)-\eD{h), 



for Model II, where 

A{x) 
Dix) 



-X 



12(1 

X 



6x(3x - 2) In X + (1 - x)(8x^ + 5x - 7) 
6a; In a: — (1 — — 5a; — 2) 



4(1 

B{x) = T^rr-^ — 7^ [(6a; - 4) Inx + (1 - x){5x - 3)] 



E(x) 



12(1 

—X 



[21na; + (l-a;)((3-a;)]. 



(7) 



(8) 



(9) 



4(1 -x)3 

Analogously, the branching ratio for b ^ sg, where the gluon is on-shell (in the sense of 
a "gluon jet"), can be written as. 



Br{b sg) 



\V*Vtt\^ Sasifi) 1 



Cl^^ifiY' Br{b^ceu), 



\Vcb\^ -Kfimc/nib) ^l{'mt/Mw,IJ') 



(10) 



where 



313063 



- 0.9135V 



0.4086 



363036J 

+0.0873r7-°-^230 _ o.0571r7^o-«994 _ o.0209r/°-^456 



(11) 



Notice the explicit yU-dependence of eq. ( p^ on ^^(yu). The fi scale of this as does not have 
to be the same as that of Cg^^{fi), but we treat them as if they are the same. From eqs. (Ef) 



and (ITOl), we form the ratio 



R = 



Br{h^sg) _ 4 as{fi) Cl"{fi) ' 



(12) 



which is independent of rric/mb. 

The scale /i denotes the renormahzation scale of the effective b ^ s'y Hamiltonian. It 
should be of order rrih, but need not be exactly equal to nib, and we shall take it to be in 



be the same as the scale at which the parameter is expanded for the QCD corrections to 
the semileptonic decay rate, eq. (||). We express as{fi) in terms of its value at /x = Mz, i.e. 



In calculating the branching ratios of 6 — * 57 and b —>■ sg in eqs. @ and (p!0D, one 
needs to know the values of the scale n and the ratio rric/mh, which are not well determined. 
However, since in our case we only want to find out the maximum value of Br{b sg) that 
is still allowed, and since the ratio R does not depend on mc/mb, we keep mc/nib at some 
fixed value ~ 1/3. Similarly, we neglect the uncertainties arising from the CKM mixing 
elements in our calculations of the b —>■ s decay branching ratios. The main uncertainty is 
therefore in the scale fi. 

We study b —>■ sg and b ^ sj numerically for different sets of parameters ^ and Mh- 
We find that a smaller value of /i gives largest Br{b sg). This is demonstrated in Fig. 1 
with ^ = 1,2 and Mh = mt = 170 GeV . In Figs. 2(a) and (b), we present the branching 
ratio oib sg decay for = 170 GeV and /i = 2.5 GeV for Models I and II, respectively. 
The hatched region to the right is ruled out by the CLEO upper bound on Br{b —>■ 57) of 
eq. (1). We notice that the CLEO limit has excluded most of the parameter space in the 
C,—Mh plane for Model II. The lower bound of Br{b 57) > 1.0 x 10~^ excludes the second 
hatched region to the left in Fig. 2(a) for Model I. We further overlay (shaded) the combined 
constraints on CKM mixing matrix (since rrit is fixed here at 170 GeV) from e parameter 
in K ^ TTTT decay, B-B mixing, and the ratios iKb/Kis] and |Kj{,/V^fe| |]18|. This constraint 



the range between 2.5 to 10 GeV as used in Ref. 0,|T7|- For simplicity we have taken it to 




(13) 



in the leading logarithmic approximation where (3q = 11 — ^ 



^Nf = 23/3 {Nf = 5). 
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is rather stringent for heavy top, and is the same for both Model I and II since the top 
coupling is common in both models. The solid, dot and dash curves in Fig. 2(a) represent 
Br{b sg) being 0.1% (< 0.1% between the two solid lines), 1% and 8%, respectively in 
Model I. For sake of illustration, however, for Model II the corresponding lines in Fig. 2(b) 
are for Br{h sg) = 0.7%, 1.0% and 1.5%, respectively. From Fig. 2 we see that, if the 
top is heavy as suggested by recent observation of CDF [^, Br{b — > sg) can at most be 



of order 1% for both Model I and II. For Model I, in fact, it would be rather difficult to go 
much beyond 0.1%. 

In Fig. 2(a), we find that the branching ratio of b ^ sg goes to zero inside the two solid 
lines due to the destructive interference between the charged Higgs and the standard model 
contributions. Br{b sg) in the region to the left of the solid line increases and reaches 
the standard model value (~ 4.5 x 10^^) when ^ = 0. Note that in the ^ — plane, for 
Model I, the interesting region allowed by the CLEO data is more or less ruled out by the 
CKM constraints, while for Model II, the variation of Br{b 57) is much faster than that 
of Br{b sg). Thus, in either case, new improvements on the experimental measurement 
of the b ^ s'-f decay will imply no significant changes on the limits of Br{b sg). 

An intriguing possibility still exists for + nih < rrit < Mw + ^6 |]2D|, allowed by all 



known constraints, including b ^ s'y and B-B mixing. When top is light, substantial charged 
Higgs contributions may in fact be called for. If one takes the heavy quark production signal 
observed by CDF [|1^ seriously, it may actually be the fourth generation t' quark. In that 
case, all loop effects are subject to GIM cancellation, and can be made ineffective. The 
crucial point, however, is that t — * bH~^ overwhelms t bW* in this domain, and can allow 
the top quark to elude past searches at hadronic colliders. For Model I, assuming that the 



heavy quark seemingly observed by CDF |jl9[| does not dominate in the loop processes, we 



find that this can happen for ^ = 1/ tan/5 ~ 2. We illustrate in Fig. 3(a) the allowed region 
for Br{b sg) for Mh < rrit = 70 GeV. Although there is still no large enhancement, the 
CKM constraint is now more forgiving, and a 3% branching ratio for b ^ sg is possible. 
This cannot be considered small when compared with the standard model expectation of 



order 10 ^ 0. For Model II, b ^ S'y provides a very stringent constraint, and in particular 



it is difficult to evade the direct search for t hH^ br^u Nevertheless, combining 
the two constraints, it is found that the region ^ = 1/ tan/3 ~ 1 is allowed for having a 
light top decaying via charged Higgs (which does not decay dominantly via rz/), especially 
when one takes into account all the possible sources of errors in making estimates. We plot 
in Fig. 3(b), with same notation as in Fig. 2(b), the expected Br{b sg) that may still be 
allowed for Model II. The solid line corresponding to Br{b sg) = 0.7% now falls outside 
of the figure. We find that the maximum value of Br{b — > sg) is about 0.9%, which is indeed 
smaller than the case for Model I. 

As stressed in ref. b ^ sg at the 10% level or higher could account for the apparent 
discrepancy on Br{b eu+X) between experiment and theory. We find that this possibility 
is quite definitely ruled out, by the combined limits of 6 ^ 57 and CKM matrix, especially if 
one takes the CDF heavy quark production signal as due to the top quark. However, in case 
the top is actually light (and CDF signal is either faked or due to new heavy quarks), b sg 
could still be at 3.5% level. Since the other possibility for suppressing Br{b ^ ev + X) by 
having b ru + X a.t 10% level or higher is also ruled out by ALEPH collaboration ||2^, the 
two Higgs doublet models cannot help alleviate the inclusive semileptonic b decay problem. 
Perhaps one would have to opt for large as{Mz) (of order 0.13) and a low /i scale {e.g. 



fi ~ WI5/2) for B decay processes, as suggested by Altarelli and Petrarca pj. This would 
imply that as > 0.3 for B decay processes. 
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FIGURES 



The ratio R as a function of the scale fj, in (a) Model I and (b) Model II. The solid and 
dash curves are for ^ = 2 and 1, respectively. 

Branching ratio oib ^ sg decay for rrit — 170 GeV and /j, — 2.5 GeV for (a) Model I and 
(b) Model II. The hatched region to the right is ruled out by the CLEO upper bound on 
Br{b — > 57) , while a second region to the left in (a) is excluded by imposing a lower bound of 
Br{b S7) > 1.0x10-^. The shaded region is forbidden by constraints on the CKM matrix. 
The solid, dot and dash curves represent Br{b — > sg) being 0.1% (< 0.1% between the two 
solid lines), 1% and 8% in (a), while in (b) they correspond to Br{b — > sg) — 0.95% , 1.0% 
and 1.3%, respectively. 

Same as Fig. 2 but with rrit — 70 GeV. 
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